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MOLECULAR PHYLOGENETIC EVIDENCE FOR THE 
PARALLEL EVOLUTION OF ROCK ECOMORPHS IN THE 
NEW ZEALAND ORB-WEAVING SPIDER 
WAITKERA WAITAKERENSIS (FAMILY ULOBORIDAE) 
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ABSTRACT. The genus Waitkera is the only New Zealand representative of the family Uloboridae and 
is known from a single species, Waitkera waitakerensis. This species is found in forests of the North 
Island, where it constructs orb-webs on understory vegetation. Rock outcrops in the Northland region 
support populations of W. waitkerensis comprised of larger individuals than those found elsewhere on the 
island, including those in surrounding forests. Parsimony analyses of DNA sequences from the mitochon¬ 
drial NADH dehydrogenase subunit ND1, using Siratoba ref ernes, another basal uloborid, as an outgroup, 
did not delineate these rock-dwelling populations as a monophyletic lineage that could be regarded as a 
distinct species. A TCS analysis leads to the same conclusion, suggesting that rock-dwelling populations 
represent independently evolved ecotypes. Northland populations of W. waitakerensis are phylogenetically 
basal; indicating that the species’ range contracted northward during the Pleistocene and recolonized the 
remainder of the North Island. 
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The monotypic genus Waitkera (Opell 
1979) is New Zealand’s only representative of 
the family Uloboridae. Waitkera waitakeren¬ 
sis (Chamberlain 1946) is comprised of small 
spiders (Fig. 1) that construct horizontal orb- 
webs in understory vegetation (Fig. 2). Fe¬ 
males have a cephalothorax-abdomen length 
of 4-5 mm and males of 3-4 mm. Adult fe¬ 
male mass averages 9 mg (Opell 1999). Evi¬ 
dence of this species can also be found in their 
peaked, triangular egg sacs (Fig. 3) that are 
deposited at the edges of their webs. In the 
Northland region, egg sacs are first produced 
in early to mid January and spiderlings begin 
to emerge from them about one month later 
(Opell unpublished observations). There are 
few published records of W. waitakerensis 
(Forster 1967; Forster & Forster 1999). How¬ 
ever, this species is sometimes the most nu¬ 
merous orb-weaver in a forest (Opell unpub¬ 
lished observations) and I have collected it 
from localities throughout the North Island 
(Fig. 4) in kauri-podocarp-hardwood, lowland 
podocarp-hardwood, and lowland hardwood 
forests. I did not observe W. waitakerensis in 
beech forests of the Huiarau or Ruahine Rang¬ 
es, the extent of my searches in this habitat. 

Pleistocene glaciation affected populations 


of New Zealand’s terrestrial arthropods (Trew- 
ick 2001) and is likely to have impacted W. 
waitakerensis . The restriction of this species 
to New Zealand’s North Island suggests that 
during the Pleistocene (New Zealand Wanga¬ 
nui series) this, and perhaps other species of 
the genus, may have been eliminated from the 
South Island. In fact, during the Pleistocene, 
it is probable that the range of W. waitaker¬ 
ensis contracted to the warm climate forests 
that persisted only in the Northland region 
(Suggate 1978; Thornton 1985). Other events 
may have contributed to the extirpation of 
these spiders. The eruption of the North Is¬ 
land’s Taupo Volcano 20,000 years ago 
(Thornton 1985) and again about 1855 years 
ago (Wilson & Walker 1985) formed deep ash 
fields that extended for hundreds of kilometers 
(Fig. 4). These eruptions have impacted pop¬ 
ulations of other species (McDowall 1996; 
Morgan-Richards et al. 2000, 2001) and 
would have been catastrophic for orb-weaving 
spiders living in the central portion of the 
North Island. 

It is possible that additional species of 
Waitkera survived these events by occupying 
refugia (Pielou 1991; Pfenninger et al. 2003; 
Trewick 2001). A likely place for these refu- 
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Figures 1-2 .—Waitkera waitakerensis : 1. Adult 
female on web; 2. Horizontal orb-webs. 

gia is the Northland region, which has been 
collected less thoroughly than most regions of 
New Zealand and in which cryptic species 
have been discovered (Gleeson et al. 1999). 
In this region I discovered a population of 
Waitkera that appeared to represent such a rel¬ 
ict species. These spiders were living in the 
cool, shaded rock crevices of the Waro Lime¬ 
stone Reserve, near Hikurangi (Fig. 5, popu¬ 
lation 4; Hawley 1981). These crevices ex¬ 
tended deep into this karstic area, as 
evidenced by the cool air coming from many 
of them. At this site, the spiders, their webs, 
and their egg sacs were conspicuously larger 
than those of forest-dwelling populations else¬ 
where on the North Island, including those in 
the vicinity of this reserve. 

The Waro rock formation is the textbook 
example of an extensive formation of soft Ol- 
igocene limestone that, after being uplifted by 
the collision the Indo-Australian and Pacific 
plates, slid westward to form the surface of 
the eastern one-third of the Northland region 
(Fig. 5; Suggate 1978; Thornton 1985). I 
found populations of the larger, rock-dwelling 
Waitkera in outcrops of this formation at three 
additional sites: Waiomio Glow Worm Cave 
near Kawakawa, Abbey Caves east of Whan- 
garei, and Waipu Caves south of Whangarei 
(Fig. 5, populations 19, 6, and 8, respectively). 



Figure 3.—Egg sacs of W. waitakerensis attached 
at the edge of an orb-web. 


1 73 * 
! 



Figure 4.-—Sites at which W. waitakerensis was 
observed by B. Opell. Numbers refer to localities 
given in Table L Heavy dashed lines denote the 
boundaries of ash deposited by the eruption of the 
Taupo Volcano about 20,000 years ago and light 
dashed lines for an eruption about 1855 years ago 
(Thornton 1985; Wilson & Walker 1985). 
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Figure 5.—Oligocene limestone distribution in 
the Northland and localities of rock-dwelling and 
forest-dwelling populations in this region. Numbers 
refer to localities given in Table 1. 


There, appears to be no mechanism for direct 
gene flow between these rock-dwelling pop¬ 
ulations. A distance of about 12 km separates 
the closest rock-dwelling sites. W. waitaker- 
ensis does not appear to be troglophilic. I did 
not find these spiders either at the entrances 
of Waipoua Cave, Abbey Caves, or Waipu 
Cave or in the first 20 m of these caves. 

Ecological and morphological evidence 
suggest that the rock-dwelling Waitkera pop¬ 
ulations may represent one or more undescri- 
bed species. However, more convincing sup¬ 
port for this hypothesis would come from 
documentation that these populations form a 
monophyletic lineage. In this case, such evi¬ 
dence can best come from a molecular com¬ 
parison of these populations. Therefore, I ex¬ 
amined the relationships among Waitkera 
populations using DNA sequences of the mi¬ 
tochondrial NADH dehydrogenase subunit 
ND1. This rapidly evolving gene has been 
used successfully to reconstruct the phyloge- 
nies and examine the population structure of 
other closely related spider species (Hedin 
1997a, 1997b; Hedin & Maddison 2001; Mas- 
ta & Maddison 2002). 

METHODS 

Morphology.—To document the morpho¬ 
logical differences among populations, I used a 


Wild M8 dissecting microscope equipped with 
an ocular reticle to measure the first femur 
lengths of ten alcohol-preserved adult female 
specimens from each of nine forest-dwelling 
populations and four rock-dwelling populations. 
I compared the mean femur lengths of these 
fourteen populations using a R y an - Ei no t-G ah ri - 
el-Welsch Multiple Range Test (alpha = 0.05) 
(Day & Quinn 1989) performed with the SAS 
statistical program (Cary, North Carolina) run 
on a personal computer. W. waitakerensis fe¬ 
males have a large, spherical median sperma- 
theca (Opell 1979). To examine differences in 
reproductive morphology relative to somatic 
features, I measured the carapace lengths and 
spermathecal widths of 24 forest-dwelling fe¬ 
males from localities throughout the North Is¬ 
land and of five females from each of the four 
rock-dwelling populations. Carapace length was 
measured under a dissecting microscope. Sper¬ 
mathecal width was measured under a com¬ 
pound microscope from a genital region that had 
been removed, cleared in clove oil, and tem¬ 
porarily mounted on a microscope slide under a 
cover slip. Voucher specimens from each of 
these populations have been deposited in col¬ 
lections of Landcare Research, Auckland and 
the Otago Museum, Dunedin, New Zealand. 

Molecular.—The molecular study included 
40 specimens from 18 populations: 3 North¬ 
land rock populations, 5 Northland forest pop¬ 
ulations, 2 populations from the north central 
region of the island, and 8 populations from 
the central and southern regions of the island 
(Table 1; Fig. 4). I used two specimens of Sir- 
atoba ref evens (Muma & Gertsch 1964) from 
the vicinity of Portal, Cochise County, Ari¬ 
zona as an outgroup for W. waitakerensis . 
This was the most basal uloborid (Coddington 
1990) for which I could obtain DNA. Voucher 
specimens from each of these populations 
have been deposited in the collections of the 
Otago Museum, Dunedin, New Zealand. 

I extracted DNA from these alcohol-pre¬ 
served specimens using a Puregene DNA iso¬ 
lation kit from Gentra Systems, Inc. and used 
PCR to amplify the double-stranded ND1 sub¬ 
unit of mitochondrial NADH dehydrogenase, 
employing the primers and thermocycler param¬ 
eters described by Hedin (1997a). For each 
specimen, I obtained the sequences of both 
DNA strands using cycle sequencing, Biosys¬ 
tem s’ Big Dye Terminator® chemistry, and a 
3100 genetic analyzer instrumentation from Ap- 
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HU forest-dwelling 



Population number 


Figure 6.—Histogram of adult female first femur 
lengths of rock-dwelling and forest-dwelling pop¬ 
ulations, arranged (from left to right) in a north to 
south order. Population numbers refer to localities 
shown in Figures 4 and 5 and described in Table 1. 
The sample size for each population is 10 individ¬ 
uals. Error bars denote ± 1 standard deviation. Let¬ 
ters within histogram bars denote the ranking of 
mean values assigned by a Ryan-Einot-Gabriel- 
Welsch Multiple Range Test (alpha = 0.05) (Day 
& Quinn 1989). 


plied Biosystems. Each electropherogram was 
edited with the EditView program and then 
aligned with its complementary sequence. Dis¬ 
crepancies were reedited and new sequencing 
reactions were run when required. Clustal V 
(Higgins et al. 1996), as implemented by DNA 
Star, was used to align sequences. 

Maximum parsimony and maximum likeli¬ 
hood analyses were implemented with Paup* 
4.0bl0 (Swofford 1998). Modeltest 3.6 (Posa¬ 
da & Crandall 1998) was run using the likeli¬ 
hood scores generated by PAUP under the 
“longfmt^yes” option of the Lscores com¬ 
mand to determine the preferred maximum 
likelihood model. The selected model and its 
values were entered into the maximum likeli¬ 
hood settings of PAUP. The TCS program 
(Clement et al, 2000) was used to perform 
Templeton, Crandall, Sing analysis (statistical 
parsimony analysis; Crandall et al. 1994; Tem¬ 
pleton 1998; Templeton et al. 1992). The 
GeoBis 2.2 program (1999-2004 David Posa¬ 
da) and its June 2004 inference key (Templeton 
et al. 1995; Posada & Templeton 2000) were 
used to perform a nested haplotype analysis 
(Templeton et al. 1987; Crandall 1994, 1996). 

RESULTS 

Morphology,—Spiders from the four rock¬ 
dwelling populations were larger than those 
from the ten forest-dwelling populations (Fig. 
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Figure 7.—Plot of adult female W. waitakerensis 
spermatheca width and carapace length. Population 
numbers refer to localities shown in Figures 4 and 
5 and described in Table 1. 


6). These size differences are also reflected in 
male and female genitalia, although there are 
no qualitative shape differences, perhaps due 
to the very simple nature of the genitalia 
(Opell 1979). As Figure 7 illustrates, rock¬ 
dwelling females have larger spermathecae 
than forest-dwellers. This may simply be a re¬ 
flection of the larger body size of rock-dwell¬ 
ers, although recent molecular studies of mil¬ 
lipedes and spiders show that genitalic 
morphology can underestimate species diver¬ 
sity (e.g., Bond & Sierwald 2002; Bond et al. 
2001). It is interesting to note that spermathe¬ 
cae of three of the four rock-dwelling popu¬ 
lations are as distinct from one another as they 
are from the forest-dwelling populations. This 
could suggest that even greater diversity is 
represented among the rock-dwellers or that 
these populations exhibit more phenotypic 
plasticity than do forest-dwelling populations. 

Molecular analysis.—ATter uncertain ter¬ 
minal DNA regions were eliminated, 412 base 
pairs were available for analysis. Three inser¬ 
tions in each W. waitakerensis sequence were 
necessary to align the DNA of this species 
with that of 5. refernes. No additional inser¬ 
tions or deletions were necessary. Sequences 
of both S. refernes specimens were identical. 
Sixteen DNA haplotypes with a mean uncor¬ 
rected (p) distance of 0.94% (range 0.24- 
1.71%) were represented by the 40 W. wai¬ 
takerensis specimens (Table 1). The distance 
between W. waitakerensis haplotypes and S. 
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Figure 8.—Parsimony analysis of 16 W. waitakerensis ND1 haplolypes run with S. referens as an 
outgroup and mapped onto localities (left) and statistical parsimony (TCS) network of these haplotypes 
(right). Parsimony analysis is the single tree from a branch and bound search. Cl = 0.9846, RI = 0.9231, 
number of informative characters = 1 16 (6 within W. waitakerensis ), 50% majority-rule consensus boot¬ 
strap values from a 2000-replicate analysis. In the TCS network open circles represent the haplotypes and 
cross-hatched circles inferred missing haplotypes. Each line connecting a haplotype represents a single 
change in one base pair. 


referens was 27.6-28.3%. GenBank accession 
numbers: Siratoba = DQ026788, HI = 
AY974175, H2 = AY974176, H3 = 

AY974177, H4 = AY974178, H5 = 

AY974179, H6 = AY974180, H7 = 

AY974181, H8 = AY974182, H9 = 

AY974183, H10 = AY974184, Hll = 

AY974185, H12 = AY974186, H13 = 

AY974187, H14 = AY974188, H15 = 

AY974189, H16 = AY974190. 

Maximum parsimony analysis (branch and 
bound search) produced a single tree (Fig. 8). 
Modeltest selected a HKY model and provid¬ 
ed the following values used in PAUP’s like¬ 
lihood settings: Ti/tv ratio = 2.6311, empirical 
base frequencies, among-site rate variation = 
0, variable sites = equal rates for all sites, 
maximum number of branch-length smooth¬ 
ing passes = 20, parameterization for clock 
model = standard, starting branch lengths for 
non-clock models = least squares method, 
Jukes-Cantor. A heuristic search produced a 
single tree (-Ln likelihood = 994.95) that was 
identical to that shown in Figure 8. An at¬ 
tempted branch and bound search was termi¬ 
nated after four days of computing. 

Parsimony analysis (Fig. 8) roots the tree at 
two haplotypes (H7 & H8) from Waipoua For¬ 
est (locality 3) in the west of the Northland. 


Relationships among the remaining Northland 
haplotypes (HI-6 & H9) are not resolved, 
providing no support for the hypothesis that 
rock-dwelling populations constitute a distinct 
lineage. Furthermore, individuals from the 
three included rocky sites share haplotype HI 
with individuals from Northland forest sites. 
Two Coromandel Peninsula haplotypes (H10- 
11) lie at the base of a clade that is sister to 
the unresolved northland haplotypes and ter¬ 
minates in five unresolved central and south¬ 
ern haplotypes (HI2-16). The only two hap¬ 
lotypes not explained by a strict cladograrn 
base-to-tip = north-to-south pattern are H9 
and HI6. H9 is represented both at locality 9, 
which is consistent with the general pattern, 
but also at locality 12, a more southern site. 
HI6 is represented at locality 9, a more north¬ 
ern locality than those occupied by the four 
other members of its unresolved, otherwise 
southern clade. 

Siratoba referens is too distantly related to 
W. waitakerensis to be included in the 95% 
confidence limits of the statistical parsimony 
network (Fig. 8). However, if this network is 
rooted at haplotype H8, its pattern is identical 
to that of the parsimony analysis. A nested 
haplotype analysis of geographical distances 
based on this network (Fig. 9) and distribution 
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Figure 9.—Haplotype network arranged in hier¬ 
archical format for nested haplotype analysis of 
geographical distances. Numbers within brackets 
identify the nested clades. Other numbers identify 
the haplotypes, shown in Fig. 8 and Table 1. 

data given in Table 1 was inconclusive. It was 
unable to determine if the current distribution 
of W. waitakerensis is explained by past frag¬ 
mentation of a more inclusive distribution, by 
continuous range expansion, by long distance 
dispersal, or by combinations of these events. 
Exact contingency tests for all subclades were 
insignificant (X 2 > 7.98, P > 0.16). The exact 
contingency test for the total cladogram was 
significant ( X 2 — 37.95, P < 0.001). For this 
clade (comprised of subclades 2-1 and 2-2) 
the inference key led through the following 
couplets: 1, 19, 20 (here, I considered sam¬ 
pling adequate to answer “yes”; however, an 
answer of “no” leads to the conclusion “in¬ 
adequate geographic sampling”), and 2, 
where the conclusion is “tip/interior status can 
not be determined—inconclusive outcome”. 
This GeoDis analysis placed the geographic 
center of clade 2-1 in the Northland, just 
north of Whangarei (between sites 4, 5, and 
7) and the geographic center of clade 2-2 
northwest of Lake Taupo (just east of site 11). 

DISCUSSION 

Phylogenetic analyses offer no support for 
the hypothesis that the rock-dwelling popula¬ 
tions of large-bodied W. waitakerensis repre¬ 
sent a monophyletic lineage that might be 
considered a species. Some of the individuals 
from each of these populations share a hap¬ 


lotype with individuals from forest popula¬ 
tions. Two unique haplotypes were included 
in one of the rock-dwelling populations, but 
these were no more divergent than other 
unique haplotypes found among Northland 
forest-dwellers. Thus, rock-dwelling popula¬ 
tions appear to represent an ecotype that has 
arisen multiple times. These populations of W. 
waitakerensis may represent lineages that 
have recently adapted to karstic areas and 
have had insufficient time to diverge geneti¬ 
cally. As such, they may represent a group of 
incipient or cryptic species. ND1 DNA may 
not have diverged rapidly enough to serve as 
an appropriate genetic marker and other genes 
might better delineate individuals from rock¬ 
dwelling populations. Adult rock-dwelling 
and forest-dwelling populations do not appear 
to be temporally isolated, although their mi¬ 
crohabitats tend to isolate them physically. 
Genitalic size differences may also limit gene 
flow between rock and forest populations, but 
the simple genitalia of W. waitakerensis (fe¬ 
males are functionally haplogyne; Opell 1979) 
may minimize the difficulty of mating by in¬ 
dividuals of dissimilar sizes. 

The larger size of rock-dwelling spiders in 
the Northland does not initially appear to con¬ 
form to intraspecific size dines that charac¬ 
terize many terrestrial arthropods with univol- 
tine life cycle. Larger individuals are typically 
found in warmer regions at the lower latitudes 
or lower elevations of a species’ range (e.g., 
Masaki 1978; Mousseau & Roff 1989; On 
1996; Schoener & Janzen 1968; Scott & Din¬ 
gle 1990). These warmer regions have longer 
growing seasons that permit individuals to at¬ 
tain larger adult sizes (Mousseau 1997), a pat¬ 
tern that appears to have an underlying genetic 
component (Huey et al. 2000). These obser¬ 
vations suggest that it may be more appropri¬ 
ate to consider rock-dwelling populations of 
W . waitakerensis as living in more thermally 
stable microhabitats rather than simply in 
cooler microhabitats. Under this scenario, spi¬ 
ders in karstic microhabitats would be pro¬ 
tected from the low and fluctuating spring 
temperatures experienced by forest-dwellers 
and, consequently, would have a longer de¬ 
velopmental period permitting them to attain 
a larger adult size than forest-dwellers. 

The population structure of W. waitaker¬ 
ensis (Fig. 8) is consistent with a contraction 
of the species’ range to the Northland region 
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during the Pleistocene glaciation. The most 
basal haplotypes (H7 & H8) are found in the 
Northland, as are seven haplotypes derived 
from them. As the climate warmed, this spe¬ 
cies recolonized the North Island by way of 
the Coromandel Peninsula (haplotypes H10- 
11). Access to more direct central and western 
corridors may have been blocked by extensive 
volcanic activity in the Auckland area (Sug- 
gate 1978; Thornton 1985). However, this cor¬ 
ridor later opened, allowing southern migra¬ 
tion of haplotype H9 and northern migration 
of haplotype H16 into the species’ type local¬ 
ity in the Waitakere Mountains north of Auck¬ 
land (locality 9). This species may have bal¬ 
looned to the Coromandel Peninsula or it may 
have expanded southward along with warm 
climate coastal forests. Colonization of the 
North Island’s central and southern regions 
appears to have been rapid and perhaps recent, 
given the occurrence of a common haplotype 
(HI2) at six widely separated sites. The spread 
of W. waitakerensis below the 38 th parallel 
may have been halted or set back by the most 
recent eruption of the Taupo Volcano (Fig. 4; 
Wilson & Walker 1985) and populations in the 
southern half of the North Island may be of 
very recent origin. 

The phylogenetic pattern of W. waitaker¬ 
ensis may be characteristic of other New Zea¬ 
land spiders that require a mild climate and 
have limited dispersal capabilities. If so, this 
study indicates that the Northland and Coro¬ 
mandel Peninsula regions, though small in 
area, contain critical elements of New Zea¬ 
land’s biodiversity. 
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